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The stability of the cosmological event horizons (CEHs) of 
a class of non-static global cosmic strings is studied against 
perturbations of gravitational waves and massless scalar field. 
It is found that the perturbations of gravitational waves al- 
ways turn the CEHs into non-scalar weak spacetime curvature 
singularities, while the ones of massless scalar field turn the 
CEHs either into non-scalar weak singularities or into scalar 
ones depending on the particular cases considered. The per- 
turbations of test massless scalar field is also studied, and it 
is found that they do not always give the correct prediction. 

98.80.Cq, 04.20Jb, 04.40.+C. 



I. INTRODUCTION 

Cosmic strings which may have been formed in the 
early Universe have been studied extensively Q], since 
the pioneering work of Kibble || . Recently, Banerjee et 
al. H and Gregory || studied non-static global strings, 
and some interested results were found. In particular, 
Gregory showed that the spacetime singularities usually 
appearing in the static case ||] can be replaced by cosmo- 
logical event horizons (CEH). This result is very impor- 
tant, as it may make the structure formation scenario of 
cosmic strings more likely, and may open a new avenue 
to the study of global strings. 

However, our recent studies || showed that these 
CEHs in general were not stable to the perturbations of 
null dust fluid, and always turned into spacetime singu- 
larities. The singularities are strong in the sense that the 
distortion of the test particles diverges logarithmically. 

In this paper, we shall study the stability of the CEHs 
against perturbations of massless scalar field and grav- 
itational waves. Specifically, the paper is organized as 
follows: in the Sec. II we consider the perturbations of 
a test massless scalar field, while in Sec. Ill and Sec. VI, 
we consider the "physical" perturbations of gravitational 
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waves and massless scalar field, respectively. The word 
"physical" here means that the back reaction of the per- 
turbations is taken into account. The paper is closed by 
Sec. V, where our main conclusions are derived. 

The main purpose of studying perturbations of test 
massless scalar field is to generalize the Helliwell- 
Konkowski (HK) conjecture about the stability of quasi- 
regular singularities Q to the stability of CEHs. As a 
matter of fact, in j(| which will be referred as Paper I, it 
was shown that the conjecture works well and gives the 
correct predictions about the stability of the CEHs, as 
far as the perturbations of null dust fluid are concerned. 

The notations used in this paper will closely follow the 
ones used in Paper I, and to avoid of repeating, some re- 
sults given there will be directly used without any further 
explanations. 



II. THE PERTURBATIONS OF TEST MASSLESS 
SCALAR FIELDS 

By requiring that the string has fixed proper width 
and that the spacetime has boost symmetry in the (t, z)- 
plane, Gregory managed to show that the spacetime for 
a U(l) global string (vortex) is given by the metric [pf 



ds 2 - e 2A ^dt 2 - dr 2 - e 2 ^ +b ^dz 2 
-C 2 {r)d0 2 . 



(2.1) 



For the cases where b(t) = ln[cosh(/3i)], +/3t, —fit, with 
(3 being a positive constant, the metric coefficients inside 
the core of a string have the asymptotic behavior 



oMr) 



0(r o - r), C(r) ~ C + O(r - rf 



(2.2) 



as r — ► Tq , where Co is a constant [cf. Eq.(3.14) in 
Ref. 5]. It was shown that in all the three cases the 
hypersurface r — tq represent a cone- like CEH |,§. 

To study the stability of these CEHs against perturba- 
tions of massless scalar field and gravitational waves, it 
is found convenient to introduce two null coordinates, u 
and v, via the relations, 



t + R 

~\7T' 



t-R 



(2.3) 
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where 

R= J e- A ^dr = ~^\n[p(r -r)}. (2.4) 

In terms of u and v, the Gregory solutions can be cast 
in the form, 



ds z = 2e- iu ^dudv 
—e 



U W [ e v (o)dz 2 +e- v ^d0 2 } 



(2.5) 



Where M (0) = y/2(3( u — v), and 



U(o) = 



In 





cosh {^-^{ u + v )^j 

^(u - v) - lnc , b(t) 

-V2/3w-lnC , b(t) 

V2j3u-]nC , b(t) 



]n[cosh((3t)], 
-Pt. 

(2.6) 



Note that in Paper I, the perturbations were considered 
in both regions r < tq and r > tq. However, as shown 
there, the conclusions obtained in these two regions are 
the same. Thus, without loss of generality, in the rest of 
the paper, we shall restrict ourselves only to the region 
r < ro- Then, the Klein-Gordon equation </>,„ g^ v = 
takes the form 



20,™ — U(o) tV <f> lit + t/(o),u0,« = 0, 



(2.7) 



where () ;X = d()/dx. 

To study the above equation, let us first consider the 
case b(t ) = +(3t. In this case, it can be shown that 
Eq.(2.7) has the general solution 



cj)(u, v) = F{u)e^ au -^/^ + G(v), [b{t) = +0t], (2.8) 

where F(u) and G(u) are arbitrary functions of their in- 
dicated arguments. To have the perturbation be finite 
initially (t = -co), we require that the two arbitrary 
functions be finite as t — > — oo and a > 0. Then, the 
trace of the energy-momentum tensor (EMT) T M1/ for the 
test massless scalar field is given by 



T 



aF(u) 



F'(u) + 



V2 



[au+P(2u-3v)]/y/2 



G'( v ) ~ ^M e (««-^»)/^ 
V2 



(2.9) 



Similarly, it can be shown that the same conclusion is 
also true for the case b(t) — — jit. 

When b(t) = ln[cosh(/3i)], Eq.(2.7) has the general so- 
lutions 



6 n e/9(«-' ! ')/v^ 



^ [(a„e^ + 2)(a„e-^" - 2) 



1/2 ; 



(2.10) 



where {b n } and {a n } are integration constants. Project- 
ing the corresponding EMT onto the PPON frame de- 
fined by Eqs.(A3) and (A4) in Paper I, we find that the 
non-vanishing components are given by 



1(0) (0) 


- T (l)(l) 




T (0)(l) 






T (2)(2) 


= T { 3)(3) 


= e 2 ? R <j> 



~<2 j2 



(2.11) 



where 
C± = 



V2/3 2 (r -r) 



-{E±e[E 2 -/3 2 (r -r)r /2 } 



n=l 



V2(3 



[ a 2 e 20R _ 4 sinh p fe PR _ 4 ]3/2 ' 

^ \al e 2 ? R - 4 sinh ftt&* - 41 3 / 2 ' 



(2.12) 



where E is a constant. From the above expressions we 
can see that, as f -> — oo, these tetrad components 
vanish, and as R — > +oo (r — > ro), the components 
T( )(Q), Tmm, and 7} )(i) become unbounded, while the 
ones T(2)(2) anc ^ ^(3)(3) remam finite. Thus, after the back 
reaction of the perturbations of the massless scalar field 
is taken into account, we would expect that the CEH is 
turned into a spacetime curvature singularity. However, 
unlike the last two cases, the nature of the singularity 
should be a non-scalar one, since now all the scalars built 
from Tfu, are finite, for example, 



T = 7\ A = — e 



2f}R, 



lif',V 



~ Const. 
Const. 



(2.13) 



as R — > +oo. To verify whether or not the above analysis 
gives the correct prediction for the stability of the CEHs, 
let us turn to consider real perturbations, that is, taking 
the back reaction of the perturbations into account. 



which diverges as r 



where a prime denotes the 



ordinary derivative with respect to their indicated argu- 
ments. Therefore, when the back reaction of the per- 
turbations is taken into account, we would expect that 
the CEH will be turned into scalar curvature singularity, 
provided that the HK conjecture continuously holds for 
CEHs §. 



III. PERTURBATIONS OF GRAVITATIONAL 
WAVES 

In Paper I, it was noted that, although the study of 
test null dust fluid and the one of real null dust fluid all 
gave the same results on the instability of the CEHs, the 
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cause of the instability was different. For the real per- 
turbations, it was caused by the non-linear interaction 
of gravitational waves, rather than what the study of the 
test particles indicated that they should be caused by the 
back reaction of perturbations of null dust fluids. Thus, 
to study the role that gravitational waves can play, we 
devote this section to perturbations of pure gravitational 
waves. These perturbations are always expected to exist, 
since at the time when the strings were formed, the tem- 
perature of the Universe was very high, and the spacetime 
was filled with gravitational and particle radiation jjj . 

To study the general perturbation of gravitational 
waves, it is found difficult. In the following wc shall study 
some particular cases. This does not lose any generality, 
since if the CEHs are stable, they should be stable against 
any kind of perturbations. Otherwise, they are not sta- 
ble. Then, from S we can easily construct the following 
solutions to the Einstein vacuum field equations, 



ds 2 = 2e~ M dudv - 



,-Ur„V r 



(e v dz z + e- y d6 A ), (3.1) 



where the metric coefficients are given by 
M = - hx[a'(u)b'(v)] - S[a(u) - b(v)] 

J_[ a{ u) + b(v)f + M c , 

V = ln[a(u) + b(v)} + 5[a{u) - b(v)} - 2mC , 
U = -la[a(u) + b(v)], 



(3.2) 



where a(u) and b(v) are arbitrary functions, and 6, Cq 
and M c are constants. The corresponding Kretschmann 
scalar is given by 



TZ — R a fj ia R a ^ lc ' 
= S 4 [12 -6 2 [a{u) 



b{v)f 



xe 



-2[S(a-b)+S 2 (a+b) 2 /4-M c ] 



(3.3) 



Choosing the null tetrad, 



in. 



-U/2 



= e- u ' 2 



ie- v ' 2 5t 



(3.4) 



we find that the non-vanishing components of the Weyl 
tensor C^\ a are given by 



5 2 b'{v) 2 e M 



{3-6[a(u) + b(v)}}}, 



1 



^2 = --C^xs [l^l x n 5 -l 



l n v m x m s ] 



S 2 e M 



a'(u)b'(v), 



*4 = 



5 2 a'(u) 2 e M 



{3 + 6[a(u) + b(v)}} . 



(3.5) 



The reason to project the Weyl tensor to the null tetrad 
is that now all the components have their direct phys- 
ical interpretations |&l9: \&o represents the transverse 
gravitational wave component along the null direction 
Ifj,, ^2 the Coulomb-like component, and ^4 the trans- 
verse gravitational wave component along the null direc- 
tion rip,. Since {rip,) defines the outgoing (ingoing) null 
geodesies p| , (^4) now represents the outgoing (in- 
going) cylindrical g ravi tational wave component. 

To use solutions (|T^) as the perturbations of gravita- 
tional waves to the Gregory solution, we have to recover 
them under certain limits. To find such limits, let us 
study the three cases bit) = ln[cosh(/3f)], +f3t, —(3t sep- 
arately. 

a) b(t) = fit: In this case, if we choose 



b{v) =Coe^ v , 



(3.6) 



and replace the null coordinate u by it', where du' — 
e~ y/ 2 ~ l3u d u, it can be shown that the solutions given by 
Eq.( |3.2[ ) reduce to the corresponding Gregory solu tion, 
as S,a(u') 0. Submitting Eq.(|J) into Eqs. (|]| and 
( |3.5[ ), we find that the Kretschmann scalar and the \1/a's 
are all Unites as t — > —00, while near the CEH where 
r — > Tq , the Kretschmann scalar TZ and the components 
\I/o and ^2 are finite, but ^4 becomes infinite. It can 
be shown that now all the fourteen scalars built from 
the Riemann tensor are finite as r — > r$ . Therefore, the 
perturbations of the gravitational waves in this case do 
not turn the CEH into a scalar singularity, although they 
do turn it into a non-scalar one. The latter can be seen 
by considering the tidal forces, represented by the tetrad 
components of the Riemann tensor in a free-falling frame 
(PPON). For example, the component i?(i)(2)(i)(2) m the 
PPON frame defined by Eqs. (A3) and (A4) in Paper I 
diverges like, 



#(i)(2)(i)(2) -»• (r-r )■ 



(3.7) 



as r 



) . Therefore, the perturbations due to the grav- 
itational waves turn the CEHs into non-scalar curvature 
singularities. However, different from the perturbations 
of null dust fluids || , now the singularity is weak in the 
sense that the distortion, which is equal to the twice in- 
tegral of the tidal forces, is finite as r — > 7yT 



'o 1 



-R(l)(2)(l)(2)^ T ^ T 



(ro-r)ln(ro-r), (3.8) 



where tq is a constant and chosen such that r — ► To as 
r -> ?v 

It should be noted that in using the PPON frame de- 
fined in Paper I to obtain the above expressions, we have 
assumed that the gravitational wave perturbations are 
weak, so the PPON frame of the perturbed solutions can 
be replaced by the one of non-perturbed solutions. This 
is the case when S, a(u') and its derivatives are all very 
small. In the following, whenever we use this frame, we 
always assume that the corresponding conditions hold. 
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b) b(t) — —0 t: In this case, to have the solutions 
given by Eq. (|3.2| ) reduce to the corresponding Gregory 
solution, as 5, b(v) — » 0, we have to choose 



a(u) — Cqc 



(3.9) 



Onc e this is done, it can be shown from Eqs^jy?) and 
(3.5) that now the CEH is also not stable and turned into 
a non-scalar singularity in a manner quite similar to that 
in the last case. In particular, the project of the Riemann 
tensor onto the PPON frame defined in Paper I diverges, 
for example, th e co mponent -R(i)(2)(i)(2) diverges exactly 
like that of Eq. ( |3.7| ) , while the twice integral of it is given 
by Eq.( |3.8| ). Thus, the non-scalar singularity is also weak. 

c) b(t) = ln[cosh/3t]: In this case, it i s ea sy to show 
that as S — > 0, the solutions given by Eq.(3.2) reduce to 



the corresponding Gregory solution, provided the func- 
tions a(u) and b(v) are chosen such that 



(3.10) 



Inserting the above expressions into Eqs.(3.3) and ( |3. 
we find that all of the fourteen scalars built from the Rie- 
mann tensor are finite both at the initial t = — oo and as 
r — ► Tq . Thus, similar to the last two cases, the perturba- 
tions of the gravitational waves do not turn the CEH into 
a scalar singularity. To see whether or not they produce 
non-scalar singularities, we can project the Riemann ten- 
sor onto the PPON frame defined in Paper I. After doing 
so, we find that some of the tetrad components indeed 
diverge, for example, the comp one nt -Rm(2)(i)(2) diverges 
exactly like that given by Eq.(3.7). 

Therefore, in all the three cases the gravitational per- 
turbations turn the CEHs into spacetime singularities, 
and the singularities are non-scalar ones, and are weak 
in the sense that although the tidal forces diverge, the 
distortion is finite. 



IV. PERTURBATIONS OF MASSLESS SCALAR 
FIELD 

To study perturbations of massless scalar field, we shall 
use a theorem given in flT^fl , which states as follows: If 
the solutions {M g ,U g ,V g } is a solution of t he Einstein 
vacuum field equations for the metric (3.1), then, the 
solution 



{M, U, V, 4>) = {M g - n g , V g , U g , XV g /V2}, (4.1) 

is a solution of the Einstein-scalar field equations G^ v — 
4>,^4 l ,v ~ 5/j^0.q0' q /2, where A is a constant, and 



Qg(U, V) = X 2 < ~Ug - ^ \2Ug <U Ug, V \ - Mg 



(4.2) 



For more details we refer readers to [[L2l . 



In order to use this theorem, the con ditio n U g:U U g:V ^ 
has to be true. However, from Eq.( |2.6| ) we can see 
that this is the case only for b(t) = ln[cosh(/3t)]. To 
overco me t his problem, we shall use the solutions given 
by Eq.( |3.2D with 6 = as the vacuum solutions for the 
cases b(t) ±/3t. It can be shown that in these two cases 
the corresponding solutions are flat, and can be brought 
to the forms that the corresponding Gregory solutions 

take by some coordinate transformations^ Once this is 

clear, we take the solutions given by Eq.(3.2) with S = 
as the vacuum solution {M g , U g ,V g } of th e Einstein field 
equations. Submitting them into Eq.( [4.l[ ), we find 

M = (1 + X 2 )M g + X 2 ln2 \a'(u)b'(v)\ 
X 2 

-yM*(«) + &(«)]. 
V = ln[a{u) + b(v)] -21nC 0) 
U = -ln[a(u) + 6(V)], 

4> = -p {\n[a{u) + b(v)\ - 2 In Co} , (4.3) 

where for the case b(t) — +/3i, the function b(v) is 
given by Eq.([3~6|), while the function a(u) is arbitrary. 
For the case b(t) = — j3t, the function a(u) is given by 
Eq^SJ}), while the function b(v) is arbitrary. For the 
case b(t) = ln[cosh(/3i)], the two func tions a(u) and b(v) 
are all fixed and give n b y Eq.( [3.1C| ). To consider the 
solutions given by Eq. ( [4.3D as perturbations of the corre- 
sponding Gregory solutions, we require that the constant 
A, the arbitrary function a(u) and its derivatives in the 
case b(t) = +f3t, and the arbitrary function b(v) and its 
derivatives in the case b(t) = —f3t, are all small. In par- 
ticular, when A, a(w) — > 0, these solutions reduce to the 
Gregory solution for b(t) = +(3t, and when A, b(v) — » 0, 
they reduce to t he Gregory solution for b(t) = —fit. 

From Eq.(4.3), we find that the corresponding physical 
quantities are given by 



T = T A A 



K 

*2 
*4 



Cl 



[a(u) + b(v)} 2 + x2 / 2 ' 
3c? 



X 2 [a{u) + b(v)} x2 + i, 

cib'jv) 
4a'(u)[a(u) + b(v)} 2+x2 / 2: 

£i 

12[a(u) + b(v)] 2 + x2 / 2, 

c\a'(u) 
Ab'{v)[a{u) +6(v)] 2 + A2 /2' 



(4.4) 



where c\ — A 2 2 a2 ^ e M =( 1 + A2 ). To study the asymptotic 
behavior of these quantities, let us consider the three 
cases separately. 

a) b(t) = jit: In this case, the function a(u) i s arb itrary 
but small, and the function b(v) is given by Eq.(3.6), from 
which we find that 



b(v), b'(v)~p(r -r)e 



0t 



(4.5) 
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Submitting the above expression into Eq.(4.4), we find 
that all these quantities are finite, except for $4 which di- 
verges like e _,3t /( r o — r ) both as t — > —00 and asr-> r^" . 
Note that the amplitude of the gravitational wave com- 
ponents \&o and \&4 is not completely fixed []. Thus, the 
divergence of ^4 does not really mean that the spacetime 
is singular. To clarify this point, let us first consider the 
fourteen scalars built from the Ricmann tensor, which are 
found finite for t — — 00 and r = ro- Therefore, in this 
case the spacetime is free of scalar curvature singularities 
both at the initial and on the CEH. To see if there ex- 
ist non-scalar singularities, let consider the tidal forces. 
Projecting the corresponding Riemann tensor onto the 
PPON frame defined in Paper I, we find that some of 
its components are diverge, f or e xample, the component 
-R(i)(2)(i)(2) diverges like Eg. (3.7), while the correspond- 



ing distortion vanishes like that of Eg. (3.8). Therefore 



in this case the perturbations of the massless scalar field 
turn the CEH into a weak and non-scalar spacetime sin- 
gularity. 

b) b(t) — —(3t: In this case, the function b(v) is ar- 
bitrary and small, and the function a(u) is given by 
Eq.( |3.8[ ), from which we find that 

-pt 



a(u), a'{u) ~ /3(tq — r)e 



00, 



as t — » —00, and 

a(u), a\u) ~ /3(r - r)e" 0t -> 0, 



(4.6) 



(4.7) 



Substituting Eg. (4.6) into Eg. (4.4) we find [1 



as r — > ? 

that all these quantities are finite as t — > —00. That is, 
the spacetime is initially free of sp acet ime si ngu larities. 
However, the combination of Eqs.(4.4) and (4.7) shows 
that the spacetime may be singular on the CEH r = r^. 
A closer study shows that this is indeed the case, and 
similar to the last case, the nature of singularity is a 
weak and non-scalar one. 

d) bit) = ln[cosh(/3t)]: In this case the two functions 
a(u) and b(v) are give n by Eq.(BTC). Combining this 
equation wit h E q.(4.4) we find that all the quantities 



given by Eq.(4.4) are finite at the initial, and diverge as 
r — > ro. That is, in this case the CEH is turned into a 
scalar curvature singularity. 



V. CONCLUSIONS 

The stability of the CEHs appearing in Gregory's non- 
static global cosmic strings [[| have been studied. It has 



in Eq.(3.4) are not completely fixed by the conditions I l\ = 
n x n\ — and n x l\ — 1. In |ld|, it was shown that in general 
they take the form Z M = B5g and = A6£, where AB = e M . 



Then, it was found that #0 = B 



(0) 



*4 



A z ~^^\ where ty^'s are independent of the choice of 



v(0) 



*2 = AB¥°\ 



and 



been shown that the gravitational wave perturbations al- 
ways turn the CEHs into non-scalar weak spacetime cur- 
vature singularities, where "weak" means that although 
the tidal forces become unbounded, the distortion re- 
mains finite as these singularities approach. It has been 
also shown that the CEHs are not stable against per- 
turbations of massless scalar fields, and are turned into 
non-scalar weak singularities for the cases bit) =±(3t and 
into scalar ones for the case b(t) = ln[cosh(/3t)]. These re- 
sults are not consistent with the ones obtained by study- 
ing the perturbations of the test massless scalar field. In 
particular, in the last case the latter predicted that the 
singularity should be a non-scalar one. Thus, to general- 
ize the HK conjecture to the study of the stability of 
CEHs, more labor is required. 
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